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ABSTRACT 

If the photospheres of solar-type stars represent the composition of circumstellar disks from which 
any planets formed, spectroscopic determinations of stellar elemental abundances offer information 
on the composition of those planets, including smaller, rocky planets. In particular, the C/0 ratio is 
proposed to be a key determinant of the composition of solids that condense from disk gas and are 
incorporated into planets. Also, planets may leave chemical signatures on the photospheres of their 
host stars by sequestering heavy elements, or by being accreted by the stars. The presence, absence, 
and composition of planets could be revealed by small differences in the relative abundances between 
stars. I critically examine these scenarios and show that (i) a model of Galactic chemical evolution 
predicts that the C/0 ratio is expected to be close to the solar value and vary little between dwarf 
stars in the solar neighborhood; (ii) spectroscopic surveys of M dwarf stars limit the occurrence of 
stars with C/0 > 1 to < 10“^; and (hi) planetesimal chemistry will be controlled by the composition 
of oxygen-rich dust inherited from the molecular cloud and processed in a dust-rich environment, 
not a gas with the stellar composition. A second generation of more reduced planetesimals could 
be produced by re-equilibration of material with dust-depleted gas. Finally, I discuss how minor 
differences in relative abundances between stars that correlate with condensation temperature can be 
explained by dust-gas segregation, perhaps in circumstellar disks, rather than planet formation. 
Subject headings: stars:abundances - planets and satellites:gaseous planets - planets and satel- 
litesiformation - methodsispectroscopic 


1. INTRODUCTION 

Ground-based Doppler radial velocity surveys and 
space-based photometric surveys have established that 
planets are very common, and perhaps ubiquitous, 
around Sun-like stars. In particular, analysis of detec¬ 
tions by the NASA Kepler mission show that planets 
with radii of 1-2 Ear th radii (Rm) are far more common 
than larger planets (|Petigura et al.l l2013t iSilburt et al.l 
[Ml). Gomparisons between estimates or upper limits 
on mass from Doppler radial velocity measurements and 
radii from Kepler also suggest that the densities of plan¬ 
ets (at least planets on short-period orbts) smaller than 
1.5M0 are consistent with a “rocky” (silicates -I- metals) 
composition, while larger planets have an a dditional low- 
mole cular weight (e.g. hydrogen) envelope (|Marcv et al.l 
l20I4fl . However, the data are insufficiently precise to 
inform as to the exact composition of these planets and 
regardless there are degeneracies in mass and radius with 
different bulk compositions. 

The photospheres of dwarf stars on the main sequence 
are broadly representative of the material from which the 
star and any circumstellar disk formed. Additional infor¬ 
mation about planet composition might therefore be in¬ 
ferred by determining elemental abundances in the host- 
star photosphere. This connection is supported by the 
observation that the abundances of refractory elements in 
the solar photosphere are, to a large extent, mirrored by 
that in primitive chondritic meteorites, the latter widely 
used as an analog of the primordial building blocks of 
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planets (e.s.. lLodder^l200^ . To a lesser extent, these 
abundances are reflected in the bulk composition of the 
planets such as Earth, but with some important excep¬ 
tions. 

The two most abundant heavy elements in the Galaxy 
and Sun (and presumably planet-forming disks) are car¬ 
bon (C) and oxygen (O). In the cool interstellar medium 
(ISM), star-forming regions, and molecular cloud cores, 
C and O are presen t primarily as the CO molecule (e.g., 
iBolatto et al.ll20I^ . Because of the unit stochiometry of 
the CO molecule, whether C or O is in excess controls 
which element is available to form other compounds and 
thus the chemistry of the gas and solids that condense or 
equilibrate with the gas. The solar ratio of these two ele¬ 
ment s is presently estimated at 0.55 ±0.12 ([Caffau et al.l 
I2011L see also Asplund et al. 2009), however this value 
is not necessarily universal. The two elements differ in 
their predicted pathways of stellar nucleosynthesis and 
C/0 has a strong positive correlatio n with metallicity 
(jTeske et al.ll2ni4lNissen et al.ll2ni4[l . 

M otivated by condensatio n scenarios for the solar neb¬ 
ula (iBarshav fc Lewis! 1197611 and obse rvations of widely 
varyi n g C/0 a i nong solar-type stars (lEdvardsson et al.] 
BOOill l. iGaido^ (j2000H proposed that C/0 controls the 
abundance of water in planetary systems. During con¬ 
densation from a hot gas with C/0 <C 1, excess oxygen 
forms refractory silicate minerals and eventually H 2 O. 
However, in systems with C/0 near a critical value 
(about 0.88), no O for water is available after silicate con¬ 
densation and planets accrete from “dry” material. At 
still higher C/0, the deficiency in O causes carbides and 
graphite to become more stable and replace silicates in 
the condensation sequence; thus “carbide planets” with 
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very un-Earth-like mineralogies, interior structures, and 
atmospheres might form around stars with C/0 > 1. 

Many analyses of the spectra of solar-type stars ob- 
tained for galactic che mical evolution (GCE) studies 
(jGustafsson et al.lll999[ i and Doppler radial velocity sur- 
yeys found that a significant fraction have C/0 > 1 
(lEcuvillon et al.l [2004 120061 : iDeleado Mena et "Si 120101 : 
iPetigura fc Marcvll201lD . A series of subsequent works 
have predicted the final bulk composition of planets 
by modeling (i) the sequential condensation of elements 
and formation of planetesimals from a cooling gas disk 
with an initial composition set by host star abundances; 
and (ii) the accretion of those planetesimals into plan¬ 
ets via integration of the dynamical equations of mo¬ 
tion. These scenarios have been applied to specific sys¬ 
tems where the photospheric abundances of elements 
have been estimated, variously predicting water-rich 
“ocean planets” and “carbon planets” lacking any water 
(iBoiid et al.l 120101: lElser et al.l 120121 : iCarter-Bond et aO 
l2012t iMoriartv et al.ll2014lb 

One question at the foundation of these works is 
whether C/0 really varies widely among Galactic disk 
stars in the solar neighborhood and ever approaches 
unity. Abundances of C and O are more uncertain than 
for many other elements because of the limited num¬ 
ber of usable absorption lines, confusion with lines of 
other elements, and non-LTE corrections. Recent studies 
have revised the C/0 of nearby solar-type stars, includ¬ 
ing those that host know n plane ts, downwards (jNissenl 
1201,11 : iTeske et'mi201,S^ . 120141 1. These newer studies 
find no cases where C/0 > 0.8, the thres hold where car - 
bide minerals are expected to first form. iFortnevI (l2012[ l 
has pointed out that the occurrence of C-rich systems 
in surveys of solar-type stars conflicts with the lack of 
known C-rich ultracool T dwar fs, as the latte r are easily 
distinguished by their spectra (jFortnevll2012h . 

A second issue is the assumption that the solids in 
planetary systems condense from a cooling gas of stellar 
composition in a circumstellar disk. Equilibrium conden¬ 
sation can explain some solids in primitive chondrites, 
widely considered analogs to the building blocks of plan¬ 
ets, but these solids constitute only a small fraction of 
those meteorites. Calcium aluminum inclusions (CAIs) 
and amoeboid olivene agregates, the most abundant re¬ 
fractory condensates in chondritic meteorite s, have vol¬ 
ume abundances < 10% and typically < 1% (lScottll20fF^ 
iHezel et al.ll2008h . The dominant constituents of chon¬ 
dritic meteorites are chondrules and fine-grained matrix, 
and there is compelling evidence that these are the prod¬ 
uct of incomplete sublimation and chemical alteration 
of pre-existing solids, e.g. older generations of solids 
and even pre-solar dust from the parent molecular cloud 
(iHuss et al.ll200.lHTrinQuier et ahlboPOUBurkhardt et al.l 
l2012li . This processing occurred under conditions that 
were very different from a gas of solar composition, i.e. 
an oxygen-rich environm ent and dust-to-gas ratios > 1 
(|Alexander fc Eb^l2012li . This environment could have 
been created by gravitational settling to the mid-plane of 
the protoplanetary disk. Because of this, the composition 
of planeteismals reflects the composition of dust grains 
in the molecular clouds cores that collapse to form stars 
and planets, plus subsequent thermal processing at mid¬ 
plane conditions. ISM dust was, in turn, the product of 
processing in the interstellar medium involving cycling 


of order 10^ times between denser, cooler clouds, where 
it accreted mantles of atoms, and the hotter intercloud 
phase, where shock-heated ions and high-energy radia¬ 
tion (U V, X-rays, cosm ic rays) sputtered atoms from the 
grains ()Tielensll2ni2al i. 

An equally intriguing possibility is that the compo¬ 
sitions of planetesimals and/or planets can be inferred 
from differential measurements of abundances in the pho¬ 
tospheres of stars. This approach is motivated by the 
detection of small differences in stellar abundances of 
elements that are more refractory (higher condensation 
temperature T/) a nd thus more likely to be incorporated 
into planets fe.g.. iMelendez et al.l 120091) . These differ¬ 
ences could arise either from the accretion of planetes¬ 
imals or planets onto the star, after the dissipation of 
disk gas, or the sequestration of solids into planets, and 
accretion of the dust-poor gas onto the star0 This leads 
to a prediction that these differences correlate with the 
presence or properties of planetary systems, and provides 
a potential “short-cut” to discovering planets as well as 
estimating their chemical composition. However, subse¬ 
quent analyses of large r samples of stars does not sup¬ 
port such a corr elation (jConzalez Hernandez et ^l20l,1l 
iLiu et ail 120141 1. and alternative explanations for these 
differences should be considered. 

In this work, I critically examine the hypothesis that 
C/0 varies significantly among neighboring stars in the 
Galactic disk (and hence the host stars of Doppler- 
detected exoplanets) from both theoretical (Sec. [2]), ob¬ 
servational (Sec. 131), and cosmochemical (Sec. H]) per¬ 
spectives. For the first I develop a GCE model to pre¬ 
dict the abundance of C and O over the history of the 
Galactic disk in the vicinity of the Sun. For the second, 
I use the spectra obtained in M dwarf surveys to place 
strict upper limits on the occurrence of C-rich single M 
dwarfs among recent large surveys. For the third, I com¬ 
bine a simple model of interestellar dust evolution with 
UV measurements of the depletion of heavy elements in 
the ISM to determine what controls the C/0 ratio and 
by how much it may vary. Finally, I develop an alterna¬ 
tive explanation that can explain the observed elemental 
relative abundance differences between solar-type stars, 
without resorting to possible signatures of planet forma¬ 
tion (Sec. [S]). 

2. EXPECTATIONS FOR C/O FROM GALACTIC 
CHEMICAL EVOLUTION 

Galactic chemical evolution (GCE) models predict 
changes in the abundances of eleme nts and isotope s 
with time in both stars and the ISM (|Prantzosl 1200811 . 
They combine yields of nucleosynthetic products in stel¬ 
lar winds and ejecta with the stellar initial mass function 
(IMF) and a prescription for star formation rate to esti¬ 
mate the production of elements and isotopes. These are 
integrated with descriptions of the flow of heavy elements 
back into the ISM and subsequent incorporation in long- 
lived low-mass stars and remnants from high-mass stars. 
Production varies with both stellar progenitor mass and 
metallicity, thus the changing metallicity of new genera¬ 
tions of stars must be tracked in a GCE. 

^ This was one explanation for the well-esta blished correla¬ 
tion between overall metallicity and giant p lanets <IGonzaledll998l : 
ISantos et aLl[2(in3 : IFischer &: Valentiif20nfi| i 
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Carbon and oxygen are synthesized in stars more 
massive than I.SMq and injected into the interstel¬ 
lar medium through winds and Type II supernovae 
SN): Type I SN are predicted to contribute < 0.3% 
Gehrz et al.lll998ll . GCE models generally predict a pos¬ 
itive trend of C/0 with metallicity (and hence time) as 
well as variation between st ellar populations with differ¬ 
ent chemical histories (e.g ., lTinslevlll98Ct iTimmes et al.l 
Il995t iCescutti et al.ll200^ . These trends are broadly ob¬ 
served within disk stars and between the thin, thick-disk 
and bulge populations, and even b etween massive and 
metal-poor irregular dwarf galaxies (iGarnett et al.lll995l: 
Cescutti et al.l 120091 iNissen et al.l I20l4 lEsteban et al.l 
20I4I1 . At least two effects are thought to contribute 
to this correlation; (i) the short main-sequence lifetime 
of massive stars, which are major contributors to a- 
elements including C and O, compared to intermediate 
mass asymptotic giant branch (AGB) progenitors, which 
contribute comparatively more to the C budge10; (ii) in¬ 
creased mass loss and C yield from younger, more metal- 
rich AGB stars at the expense of O. Mass-loss rates de¬ 
pends on metallicity through its effect on the opacity 
of outer stellar atmospheres and the impact on their 
structure. The implications of the abundance depen¬ 
dence of G/0 yields for the chemical evol ution of galax¬ 
ies has been previousl y pointed out (e.g., iMaeded 119921 : 
iFraver fc Brownl[l997tl . 

The G/0 of bulge and disk populations diverge due 
to a combination of differences in the stellar IMF, wind- 
driven mass loss in the former and the accretion of pri¬ 
mordial metal-poor gas in the former. The flatter IMF 
of bulge stars implies a higher relative production of O 
from more abundant massive stars. Addition (loss) of 
mass from a star-forming system will increase (decrease) 
the amount of star formation required to arrive at a given 
metal abundance and therefore the G/0 at that abun¬ 
dance. 

One feature of the statistics of stellar metallicities of 
the Galactic disk in the neighborhood of the Sun is the 
relatively narrow dispersion and paucity of metal-poor 
stars. This is often called the “G dwarf problem” and 
translates into a flat age-metallicity relation for the most 
of the history of the Galactic disk, at least in the solar 
neighborhood. It is actually not a problem for GGE mod¬ 
els if addition of metal-poor gas to the disk is admitted 
to ba lance stellar nucleosynthesis (e.g., iHolmberg et al.l 
12001 . If G/0 is strongly correlated with metallicity 
then a flat age-metallici ty relation might me an that G/0 
evolves very little le.g.. iWheeler et ahlllOSOH . 

To re-visit the question of the evolution of G/0 in the 
solar neighborhood, the abundances of the two elements 
were calculated on an isotope-by-isotope basis using the 
model described in detail in the Appendix. The model 
simulates the production of these isotopes and release 
into the ISM in SN explosions, winds from massive stars 
and AGB stars, sequestration of isotopes into long-lived 
low-mass stars, and addition of metal-poor gas by infall 
onto the Galactic disk. The ISM is described by two com¬ 
ponents; an inter-cloud medium which produces molecu- 

® While O yields are relatively insens itive to model parameter s 
and are consistent from model to model HWooslev &: Weaverlll995h , 
the C yield of intermediate-mass s tars depend sensitive ly on the 
amount of convective “dredge-up” IIHenzini &: Volilll981ll . 


lar clouds but no stars, and a star-forming giant molecu¬ 
lar cloud component. Yields for intermediate- and mas¬ 
sive stars are taken from a variety of recent sources. Best- 
fit parameters describing the timescale of exponentially- 
declining gas infall on the disk, the initial metallicity of 
the infalling gas, the age of the stellar population at the 
solar galactocentric radius, the index of the power-law 
describing the IMF for massive stars, and the power-law 
relationship between star formation rate and gas surface 
density were found by a Monte-Garlo Markov Ghain anal¬ 
ysis: The observational constraints were the age of the 
Sun, the present mass surface density of stars, stellar re- 
manants, and gas at the solar galactocentric radius, solar 
metallicity (here taken to be G-l-0), the current metal- 
licity of the ISM, and its intrinsic standard deviation. 

The predict ion evolution of G/O is compared to the so¬ 
lar value from iGaffau et al.l (|2nilh in Fig. [T] After initial 
transients that die away after a few hundred Myr, there 
is only moderate evolution in G/0 during most of the his¬ 
tory of the Galactic disk. The G/0 ratio rises from 0.5 
to about 0.65 over the first 6 Gyr in response to increas¬ 
ing input from stars of lower mass and higher metallicity, 
as described above, and declines slightly thereafter. The 
predicted value at the time of the Sun ’s formation is 0.64, 
withinin Icr of the lGaffau et al.l (1201 Ifl value. The formal 
error in the predicted value, based on the standard devia¬ 
tion of the MGMG chain after removal of the “burn in”, is 
only 0.003, but this value does not reflect the dominant 
source of uncertainty in these calculations - the nucle- 
osynthetic yields. Although the model is unable to ex¬ 
actly reproduce the solar values of G and O abundances, 
correctly re producing these using GGE been always been 
challenging. iMaeden (|1992h also found that satisfactorily 
reproducing a “flat” C/0 places constraints on the rem¬ 
nant or “cut-off” mass at the center of a SN progenitor 
which is not injected into the ISM. 

Figure [5] plots the predicted C/0 evolution vs. pre¬ 
dicted [0/H], again co mpared to the Sun. Also plotted 
are the values found bv INissen et al.l (I20I4I1 for Galactic 
disk stars (filled points) and bulge or halo stars (open 
points) and converted to absolute values assuming a so¬ 
lar C/0 of 0.55. The predicted trend is perfectly consis¬ 
tent with the disk values if an offset of about 16% with 
the predicted vs. measured Solar value is artificially re¬ 
moved. The absence of data for disk stars at low [0/H] 
is a consequence of the age-metallicity relationship: com¬ 
paratively little time elapsed (< 1 Gyr according to the 
model) and few stars formed in this interval. 

The predicted constancy of C/0 contrasts with pre¬ 
vious predictions for an incre ase in C/0 wi t h met allic- 
ity/time. For e xample, both ICescutti et al.l (|2009ll and 
iMattssonl (|201(Tl l predict a positive slope of ~ 1 dex/dex 
in [C/0] vs. [0/H], and thus a C/0 of about 0.06 at 
[0/H] = -1. Since the model and the best-fit parame¬ 
ter values presented here are similar to those of previous 
works the most likely cause of this difference is in the par¬ 
ticular nucleosynthetic yields that were used. T he yields 
from AGB stars have been significantly revised ([Karakasl 
12010(1 but in the model these are minority contributors 
to the C and O budgets. Indeed, if there contribution 
is completely removed, C/0 still changes l ittle and the 
predic ted solar C/0 is slightly closer to the lGaffau et al.l 
(|2011ll value. Instead, the predictions of this model prob¬ 
ably stand out from previous results because of the dom- 
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Figure 1. Evolution of C/O at the solar galactocentric radius 
according to the GCE model described in the text and Append ix. 
The solar value and its uncertainty from ICaffau et al.l 11201111 is 
plotted. 

inant contribution of massive stars to both C and O. 
Since these stars have very short lives compared to the 
chemical evolution of the disk, there is no effect from 
delayed contribution, and the C/O rapidly approaches a 
steady-state value. 

Previous works have pointed to favorable comparisons 
of model predictions with observations, but the latter 
are combined data on bulge, halo, and thick- and thin- 
disk populations, therefore assuming a common origin for 
these populations that may not exist . For example, the 
[C/01 of thi n disk stars pre s ented by iBensbv fc Feltzin^ 
!|2008ll and iCescutti et al.l (1200911 show no correlation 
with [O/H] and it is only when bulge and thick disk stars 
are added that such a correlation appears. The apparent 
consistency betwee n the predicted and observed trends 
of [C/O] vs. [0/H] (jCescutti et alJl2?)f)^lMattssonll2010ll 
belies the fact that the models were tuned for the (thin) 
disk population in the solar neighborhood and should not 
be compared with other stellar populations. 

It is indisputable that metal-poor bulge/halo stars 
have low C/O compared to the sol ar neighborhood 
(jFabbian et al.l[20?)^lNissen et al.ll2ni4[ l. The model pre¬ 
sented here cannot, and was not designed to, explain 
these populations, but their abundances might be a re¬ 
sult of a flatter IMF, early loss of gas and/or truncation of 
star formation. But bulge/halo stars are rare in the solar 
neighborhood and because their lines are very weak they 
are usually avoided by planet searches using the Doppler 
radial velocity method. 

C/O may also vary at some level within the solar galac¬ 
tocentric annulus because the sources of nuclides (stars 
and stellar clusters) are discrete and mixing by rotational 
shear and epicyclic dynamics is not completely efficient. 
However, GCE is occurring on a timescale much longer 
(few Gyr) than the rotational time of the Galaxy at 8 kpc 
(^^250 Myr) which governs the rate of mixing, and thus 
this limits the magnitude of such heterogeneities. 

3. M DWARFS AS PROBES OF HIGH G/O 

In solar-type stars, the absorption lines of C and O 
are weak and C/O has little effect on the overall spec¬ 


Figure 2. Evolution of C/O at the solar galactocentric radius 
according to the GCE model described in the text and Appendix, 
plotted vs. [O/H]. The lCaffau et aP J2Q1U) solar value is plotted as 
well as values for Galacti c disk stars (filled p oints) and bulge/halo 
stars (open points) from fNissen et ^ (|20l4 ) and tied to the solar 
value. 


trum. iFortnevI (|2012t1 pointed out that cooler stars, i.e. 
brown dwarfs, offer more obvious constraints on the oc¬ 
currence of high C/O systems. The spectra of M dwarfs, 
with effective temperatures Tes < 3900K, also include 
absorption bands of 0-containing molecules, e.g. TiO, 
VO, and Ca OH, which are sensi tive to the available O 
abundance (jSchmidt et al.l 12009( 1. and hence C/O. At 
C/O > 1, TiO is absent and the Schwan bands of C 2 , 
CN, and CH should appear. Such spectra are character¬ 
istic of carbon stars, evolved stars in which dredge-up of 
carbon-righ interiors has occurred (AGB stars). If the 
initial mass function of star formation is chemically in¬ 
variant, a limit on the occurrence of C-rich M dwarfs is 
also a limit on C-rich G dwarfs. 

Dwarf carbon (dC) star have been identified 
(|Dahn et al.l 11977(1 but these probably accreted carbon- 
rich gas from a present or form er evolved companion 


(|Behara et al.[[20ir)(: IC reed 


uncommon ( Downes et al.[ 


2015) . Overall, dC stars are 


2004 ICreed (201,31 1. Carbon 


enhanced metal-poor (CEMP) stars in the Galactic halo 
appear are also rare overall but more prevalent among 
(initially) metal-poor syst ems because less a ccreted C is 
required to increase C/O (|Beers et al.[[1992ll . Any tally 
of intrinsically carbon-rich dwarfs must remove these in¬ 
terlopers or be considered an upper limit. 

I considered three spectral indices which are 
continuum- normalized measu res of the emission in spe¬ 
cific bands ((Reid et al.[[1995|) : CaH, which is the mean 
of the CaH2 (6814-6846A) and CaH3 (6960-6990A) in¬ 
dices; Ti05 (7126-7135A); and CaOH (6230-6240A). 


Values of these temperature- and gravity-sensitive in¬ 
dices are highly and positively correlated among solar- 
metallicity M dwarfs. Carbon-rich stars can be identi¬ 
fied by relatively weak bands of oxygen-containing TiO 
and CaOH, i.e. high Ti05 or CaOH indices, for a given 
strength of the non-oxygen-containing CaH. However, in¬ 
dices of metal-poor subd warfs or “extrem e ” subdwarfs 
exhibit similar behavi or ((Jao et al.[ [2008( : (Woolf et al.[ 
[2009( : [Raiourohit et al.[[^14D . 
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Figure 3. Top: Normalized PHOENIX BT-Settl synthetic spec- 
tra of M dwarfs with Teff = 3700 K, log 5 = 5, and solar abundances 
UCaffau et aLlf^lll . red line), metal-poor subdwarf ([Fe/H] = -1, 
[a/Fe] = 0, green line), extreme subdwarf ([Fe/H] = -2, [a/Fe] = 
-1-6.4, \ogg = 5.5, blue line), and solar-metallicity star except some 
O as C such that C/0=1 (red line). Bottom: ratios of the three 
non solar-metallicity spectra described above to the solar metallic- 
ity spectrum. 

To guide discrimination between C-rich M dwarfs and 
metal-poor stars based on indices, synthetic spectra were 
generat ed using the BT-Se ttl version of the PHOENIX 
models (lAllard et al.ll20lill . Four cases were considered: 
a solar-metallicity star, a star with [M/H] = 0 but C 
and O adjusted so that C/0=1, a metal-poor subdwarf 
(sd) with [M/H] = -1, and an extreme subdwarf (esd) 
with [Fe/H] = -2. The last two cases also had an alpha- 
element enhancement [a/Fe] = -1-0.4. Effective tempera¬ 
tures over the entire M dwarf and late-K dwarf range in 
steps of 100 K were considered. All cases had logg = 5, 
except for the esd which had log(; = 5.5. The visible 
and far-red portions of the spectra are plotted and com¬ 
pared in Fig. |3l The bands used to compute the CaH, 
Ti05, and CaOH indices are indicated. The current BT- 
SETTL line lists do not actually include CaOH, but the 
bandpass also includes the TiO 7 line which is also sen¬ 
sitive to C/0. 

As expected, in the C/0=1 case (black lines in Fig. [3]) 
the TiO bands at around 6235 and 6700A are dramati¬ 
cally weakened relative to the solar case (red line) while 
those of CaH are essentially unchanged. This is also true 
for the esd case (but not the subdwarf case). The differ¬ 
ential response of Ti05 and CaH is th e basis for the “C” 
parameter developed bv iLeoine et al.l (|2007f) . However, 
the C/ 0=1 and esd cases differ markedly at 4700-5300A 
and below 4500A because Fe lines are very weak in the 
esd case. There are also difference at A > 8OOOA (Fig. 
El). Absorption by H 2 O over a broad wavelength range 
centered at 1.9/im is also weaker in the C/0=1 case, but 
this would be difficult to ascertain from the ground. In¬ 
dex values were calculated for the four cases at each Teff. 

I searched for carbon-rich M dwarfs in three spectro¬ 
scopic catalogs: the CONCH-S HELL spectroscopic cat¬ 
alog of nearby, bright M dwarfs (IGaidos et al.ll20lll , the 
Sloan Digital S ky Survey ( S PSS) catalog of M dwarfs 
constructed by iWest et 5] (1201 111 , an d the LAMOST 
spectroscopic catalog of M dwarfs in lYi et al.l (I20I4D . 
Stars in the CONCH-SHELL catalog were selected based 
on brightness (J < 9), parallaxes or proper motions con- 



Figure 4. CaH vs. Ti05 indices for 2583 stars from the CONCH- 
SHELL catalog of IGaidos et al.l <201411 . some with repeated obser¬ 
vations. The black dashed curve is a quadratic fit to the locus. 
Eighteen stars with anomalously weak Ti05 bands (large indices) 
for their CaH band strength are marked as magenta points. The 
circled point (Wolf 1130 or GJ 781) has the most deviatory Ti05 
index. The isolated point illustrates the median errors. 

sistent with main sequence status, and colors. Although 
some sets of stars in the catalog were selected based on 
optical and infrared colors consistent with “normal” M 
dwarfs, this criterion was relaxed for another set, thus 
admitting any C-rich dwarfs with peculiar colors. The 
selection for “red” {V — J > 2.7) stars is not relaxed, and 
this inevitably eliminates some metal-poor early M-type 
dwarfs, although not necessarily C-rich examples. Fig- 
urejUplots the Ti05 vs. CaH indices for CONCHSHELL 
stars. The red line is the best-fit quadratic to the locus 
and the black points and lines are the values computed 
from PHOENIX BT-SETTL spectra. 

I identified 18 CONCH-SHELL stars with a Ti05 band 
that is significantly (> 3cr) weaker (larger index) com¬ 
pared to a best-fit second-order fit of Ti05 as a function 
of CaH. When calculating the significance of a devia¬ 
tion, an intrinsic dispersion of 0.032 in Ti05 was added 
in quadrature to the formal measurement errors. None 
of these 18 stars had indices that deviate to the extent 
predicted for a C/0=1 star. The most deviant star is 
PM I20050-b5426/GJ 781/Wolf 1130, with a Ti05 in¬ 
dex that is 0.17 above the best-fit locus value. Thi s was 
previ ously identified as an active Ml.5 subdwarf (iGizisI 
II 997 I) with a me tallicity based on an inf rared spectrum 
of —0 . 64 ± 0 .17 (|Roias-Avala et ani2012D . Interestingly, 
iCizisI (I1998D identified Wolf 1130 as a single-lined spec¬ 
troscopic binary [P rs 0.5 d) and proposed that the un¬ 
seen companion was a O.3M0 helium white dwarf. A 
Hubble Space Telescope Imaging Spectrograph (STIS) 
spectrum of the star (Fig. EJ exhibits some TiO ab¬ 
sorption and no diagnostic carbon-star features. There 
is also no indication of a WD companion, although this 
would be consistent with an advanced age and hence low 
UV luminosity of any such object. Spectra of the other 
candidates indicate they are metal-poor stars, or have 
systematic errors. 

I examined the molecular indices for two much larger 
samples of M dwarfs from the SDSS and LAMOST 
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Figure 5. HST STIS spectrum of Wolf 1130/GJ 781, the star with 
the comparatively weakest Ti05 band in the catalog (circled point 
in Fig. 11. Some atomic lines and molecular bands are labeled; the 
red dashed lines mark the locations of some CN and C 2 Schwan 
bands observed in carbon stars, but not in this star. 

(jWest et al.l 120111 : lYi et al.l I2014D . Figure |6] shows the 
Ti05 and CaOH indices vs. CaH index for 70,841 M 
dwarfs with spect ra in Data Releas e 7 of the Sloan Dig¬ 
ital Sky Survey (iWest et al.l 1201 111 . The grey scale is 
linearly related to the density of stars in the Ti05- or 
CaOH-vs. CaH plots. The black dashed line is the best- 
ht to the stellar locus, and the red, aquamarine, dashed 
purple, and black solid lines are the predictions from the 
BT-SETTL models for the solar metallicity, subdwarf, 
extreme subdwarf, and C/0=1 cases, respectively. Stars 
for which either index value is significantly (> Scr) higher 
than the best-fit polynomial to the stellar locus, and 
within Icr of or above the predicted C/0=1 locus and are 
plotted as open points. The 30 stars where both indices 
satisfy these criteria are plotted as filled points. Inspec¬ 
tion of the SDSS spectra of these 30 stars found that all 
are consistent with template spectra of solar-metallicity 
or metal-poor stars and do not have the featuers expected 
of C/0=1 stellar atmospheres. 

Figure [7] is the analogous set of plots for 67,082 can¬ 
didate M dwarfs from a pilot survey of the Large Area 
Multi-Object F iber Spectroscopic Telescope (LAMOST, 
lYi et al.l I2014D . The intrinsic scatter in index values, 
after correcting for measurement errors, is much larger 
for this sample, perhaps due to systematic errors. Only 
two stars have both Ti05 and CaOH values significantly 
above the locus, and within Icr of the predicted C/0 line. 
However, the index values of these stars are all ^ 1 and 
probably spurious. The paucity of high C/0 candidates 
in the LAMOST survey compared to the SDSS DR7 sam¬ 
ple is undoubtedly due to the larger uncertainties in the 
indices. 

Based on the CONCH-SHELL sample alone, high C/0 
^ 1 stars constitute less than 1.2 x 10“^ of M dwarfs 
with 95% confidence, a s tricte r constraint than the limits 
placed bv INissen et al.l (1201411 based on spectroscopy of 
66 solar-type stars. The null result from the DR7 sample 
places an an even more stringent upper limit of 6 x 10“"^ 
at 99% confiden ce. It is possible that manual screen¬ 
ing performed bv IWest et al.l (1201111 to remove corrupted 
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Figure 6. CaH, Ti05, and CaOH indic es of M dwarfs in Data 
Release 7 of the SDSS llWest et al.lI201H) . The density of stars 
is shown as a grey scale. The dashed black curves are the best 
fits to the loci. The red, blue, purple dashed, and black solid 
curves are the predictions of the PHOENIX BT-SETTL model for 
solar-metallicity, subdwarf, extreme subdwarf, and C/0=1 cases, 
respectively. 

spectra could have excluded carbon-rich stars, but at the 
resolution and signal-to-noise of the spectra, the differ¬ 
ences that distinguish C/0=1 spectra from the solar case 
would not make them overwhelmingly abberant. 

4. FROM DUST TO PLANETESIMALS 

Theoretical studies linking the abundances in the pho¬ 
tospheres of host stars to the composition of rocky plan¬ 
ets generally assume that the building blocks of planetes- 
imals condensed directly from a cooling gaseous disk of 
the same composition, but there is accumulating mete- 
oritic evidence that primitive solids in the Solar System 
were by and large not produced by de nova condensa¬ 
tion from a gas of solar composition. Instead, this ma¬ 
terial was the outcome of partial thermal processing of 
pre-existing solids under high dust/gas ratios and oxy¬ 
gen abundance (fugacity, / 02 ) brought about by growth 
and settling of grains to the nebular mid-plane and/or 
inward mi gration of water-rich ( and hence 0-rich) plan- 
etesimals (IGrossman et alJ[2008[ l. This evidence includes 
(i) the survival of pre-solar grains and isotopic anomalies 
inherited from the mo lecular cloud and older generations 
of stars (|Davi iEnni); ( ii) the gradual variation in the 
abundance of elements with condensation temperature in 
primitive meteorites, which contradicts the sharp cutoff 
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Figure 7. Same as Figure 6, but for candidate M dwarfs from the 
LAMOST sample. Several stars with anomalous index values lie 
above the plots. 

predicted by equilibrium condensation and which canno t 
be explained by radial transport in the disk (lCieslall2008ll : 
(iii) the retention of volatile sodium and sulfur and lack 
of isotopic mass fraction ation for potassium and silicon in 
chondrules (jScottll20(T^ : and (iv) molyb denum and tung¬ 
sten depletions in refractory inclusions (iFeelev fc Palmd 
Il985f) . The equilibrium oxygen fugacities of minerals in 
many meteorites also suggest high ambient / 02 , but some 
could also be explain ed by alteration on the parent bod¬ 
ies of the meteorites (|Grossman et al.ir2008ll . 

Dust growth and settling to the midplane of a disk is 
predicted to occur in 10 ^- 10 ^ yr, depending on the inten¬ 
sity of turbu lence in the disk (|Nomura fc Nakagawal2005 
lCieslal[200^ . Dust settling, as well as growth, might be 
observed via it s effect on the spe ctral energy distribu¬ 
tion of a disk (| Tanaka et al.l[2005ll . Tentative evidence 
for significant dust settling and depletion from the upper 
layers of the disks o f T Tauri stars (a ges of ~ 10® yr) 


has been presented (iFurlan et al.l 12006 

), but unambigu 

ous detection of settling is challenging ( 

Murakawall2014ll 


If settling occurs faster than the viscou s accretion time 
of a disk (10® vr. lHartmann et'HIllOQSll then grains will 
experience high temperatures only in a dust-rich envi¬ 
ronment. 

If the precursor material of planetesimals is dust, 
rather than gas, then to a large extent the composition 
of interstellar grains governs the composition of rocky 
exoplanets. Interstellar dust begins its existence as con¬ 


densates in the cooling envelopes and winds of AGB and 
red giant branch stars but these grains are subsequently 
and completely altered by many cycles of erosion and 
formation of mantles in the ISM. Erosion takes place in 
the lower-density, higher-temperature inter-cloud phase 
of the ISM, principally by sputtering by ions heated by 
the passage of supernova shocks as well as UV pho¬ 
tons. Gondensation takes place onto surviving grains 
that are incorporated into the denser, cooler cloud phase 
of the ISM. The cycling time between these two phases 
(~ 3 X 10^ yr, set by the cloud lifetime) is short compared 
to the mean time since formation in a circumstellar wind 
(~ 3 X 10® yr), hence the bulk elemental composition of 
ISM grains is s et by the bal ance between condensation 
and sputtering (lDrainell2003ll . 

The bulk composition of dust can be inferred by mea¬ 
suring the depletion of elements from the gas phase, i.e. 
using the strength of UV absorption lines along different 
lines of sight through th e ISM to some suitable back¬ 
ground source. Following iTielensI (II998I1 . the equations 
of motion of the depletions S in the cloud (c) and inter¬ 
cloud (i) medium of any particular element can be writ¬ 
ten as: 

^ = -k2 {5c - Si) + k4{l- 5c ), ( 1 ) 

at 

and 

| = ( 2 ) 

respectively. In these equations, ki, k 2 , k^, and ^4 
are the rate of mixing from the intercloud to the cloud 
medium, the rate of mixing from the cloud to the inter¬ 
cloud medium, the rate of grain destruction in the inter¬ 
cloud medium, and the rate of grain growth from molec¬ 
ular cloud gas, respectively. The nucleosynthetic rate of 
production of an element is much slower than any of the 
rates of the formation of clouds, dissipation of clouds, 
destruction by sputtering in the inter-cloud phase, and 
growth condensation in the cloud phase, respectively, and 
is ignored here. The steady state solutions are: 

^ l + fci/fca 
" l + ki/k^ + k^/ki^ 

and 

I + ki/kz + k2/k4, 

Thus, the steady-state abundances are governed by only 
two parameters, the amount of growth in the clouds 
^ 4 /^ 2 , and the amount of erosion between clouds: k^/ki. 

I estimated parame ter values for s ome elements using 
the data compiled bv IJenkini (I2009I1 . I set 5c and 5i to 
1 — [X/H]i (maximum depletion) and 1 — [X/H]q (min¬ 
imum depletion), respectively, and solved for fca/fci and 
^ 4 /^ 2 - Figure E] plot s these two py ameters. Three el¬ 
ements with data in iJenkinsI (120091 1 are not shown: to 
explain the abundances of P, Cl, and Zn in the context 
of this model requires negative destruction in the inter¬ 
cloud medium. This may be an artifact of photoioniza¬ 
tion since the abundance of these elements is estimated 
from their singly ionized forms. 

The comparatively low condensation rates of Kr, C, 
and O reflect the volatility of these elements. The be¬ 
havior of O deviates strongly from that of the refractory 
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elements, perhaps because the primary carrier of O is 
not silicates but a much more volatile substance such 
as water ice. Lik ewise C does n ot behave as a refractory 
such as graphite (jTielensll2012bjl . This sugg ests that C i s 
present as relatively volatile organic matter (j,Tonesll2nn9[l . 
Kr is not expected to condense and the non-zero value of 
ki/k 2 may be a consequence of measure ment errors or de¬ 
partu res from solar relative abundances (iCartledge et al.l 
1200 ^ . 

The differences in the parameters for C and O man¬ 
ifest themselves as a modified C/0 ratio in interstellar 
dust with respect to the total (gas-|-dust) abundance. In 
the case of th e nominal rates inferred from the data of 
iJenkinsI (|2009ll , the degree of depletion from cold molecu¬ 
lar cloud gas shows that the C/0 of the dust is 0.92 times 
that of the bulk ISM, i.e. slightly more oxygen-rich than 
the current bulk ISM value. Although O is depleted more 
rapidly from interstellar grains, it is also accreted more 
rapidly. One caveat of this estimate for the dust C/0 is 
that the gas-phase depletion of C is more uncertain than 
other elements because there are few suitable absorption 
lines (|Sofia fc Parvathill2009ll . 

The C/0 could also vary bet ween locations as a result 
of varying k^/ki and/or ^ 4 /^ 2 . iTielensI (|I998ll estimated 
the residence time 1/^2 in the warm inter-cloud phase 
of the ISM as 3 X 10® yr assuming that it is set by the 
timescale for shocking by SN and subsequent cooling and 
collapse into clouds. This was consistent with an ISM 
mass fraction in the warm phase of ~10%. However, 
the mass in the warm phase is prob ably comparab le to 
the dense molecular H 2 (cloud) phase IPraind Ornll and 
thus the residence time in the inter -cloud phase is similar 
to the cloud lifetime 3 x 10^ yr (iMnrraH 1201 III . Over 
this time, dust grains may experience ~ 10 SN shocks 
before becoming incorporated into molecular clouds. In 
contrast, ki/k 2 would be expected to vary only to the 
extent that the lifetimes of molecular clouds varies. 

I estimated the sensitivity to variations in these rates 
by multiplying each of the parameter ratios by varying 
factors e [0,3]. Figure [S] plots the predicted variation 
of C/0 with contours intervals of 10%. The unadjusted 
parameter values, which predict [C/0]dust = 0.92 (heavy 
contour), are at unit value abcissa and ordinate (circle). 
Variation in the efficiency of grain destruction in the in¬ 
tercloud medium {k^/ki) have more effect on dust C/0 
than variation in the efficiency of grain growth in molec¬ 
ular clouds {ki/k 2 )- This difference is a consequence of 
the larger dispersion (a ratio of 30) in the removal of O 
vs. C during grain destruction compared to the disper¬ 
sion in the incorporation of O vs. C during grain growth 
(a ratio of ^ 3, Fig. |8]). The very volatile behavior of O 
compared to the other elements is presumably because 
some of it is incorporated as water ice mantles around 
dust grains. 

This model predicts that large (factors of two) varia¬ 
tion in dust lifetime produced by different shock and UV 
conditions in the intercloud medium will produce mod¬ 
est (~20%) variation in bulk dust C/0, with dust in the 
vicinity of massive SN progenitors more carbon-rich. By 
the arguments presented above, variation in dust C/0 
could generate diversity in the composition of planeteis- 
mals, but this will be limited by the extent that water 
ice is retained on dust grains (see below). 



Figure 8. Inferred rates of growth (in molecular clouds) and 
depletion (in the inter-cloud medium) of ele ments in inter stellar 
grains based on the observations analyzed in IJenkini II2009I 1. The 
rates are normalized by the residence time of grains in the cloud 
and inter-cloud media, respectively. 



Figure 9. Predicted variation of C/O in interstellar dust with 
rates of dust accretion and destruc tion relative to canonical pa¬ 
rameter values, using the model of I'fieleni II1998I1 and canonical 
values set by the observed depletions in the ISM. Contours are in¬ 
tervals of 0.1, increasing upwards, and the canonical rates yield a 
C/O of 0.92. Rates for each element are assumed to vary by a 
uniform multiplicative factor. 

5. ELEMENTAL ABUNDANCE VARIATION AMONG 
(PLANET-HOSTING) STARS 

A pronounced correlation between the metallicity of 
the host star and the occurrence of giant planets was 
discovered early in the history of exoplanet science 
(lGonzalezlll998fi . Increasingly precise measurements of 
the relative abundances of elements in stellar photo¬ 
spheres have permitted more nuanced i nvestigations of 
corre la tions with plane t occur rence fe.g.. lRobinson et al.l 
I2006D . IMelendez et al.l (I2009D found that the solar com¬ 
position was anomalous compared to 9 out of a sam¬ 
ple of II solar “twins” and that the Sun is deficient 
(< 0.1 dex) in refractory elements (condensation temper¬ 
ature Tc > 900 K) compared to more volatile elements. 
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They proposed that this was a signature of the seques¬ 
tration of t hose elements in rock y planets and cores of 
gas giants. iRamirez et al.l (1201011 also found a statisti¬ 
cal correlation between planets and a deficit of refrac¬ 
tory elements among larger samples of solar-type stars 
in Doppler surveys. If differences in relative photosphere 
abundances are a signature of rocky planets or giant 
planet cores, this would provide a short-cut to planet 
detection and even a m eans to probe the composit ion of 
the planets themselve s (iDelga do Mena et al.ll2 0inH 

On the other hand, iGonzalez Hernandez et ahl (1201,11 1 
find no clear correlation between the slope of the rel¬ 
ative abundance vs. Tc and the presence or absence 
of planets. One factor that may contribute to these 
contrasting results is the removal of galactic chemical 
evolution trend s, e.g. decreasing [g/Fe] w i th inc reasing 
metallicity, by IGonzalez Hernandez et ahl (|2ni,l[l . Mo- 
roever, lOnehag et al.l (1201411 found no trend with Tc 
in the relative abundances of members of the solar-age, 
solar-metallicity cluster M67. Observations of physical 
pairs of stars with and without (detected) planets pro¬ 
vide the clearest test of the planet-abundance relation, 
as both stars should have formed from the same molec¬ 
ular cloud. But the available obs ervations do unambigu¬ 
ously support s uch a correlation (iTucci Maia et al.ll20l4 
iLiu et ai.N2014ll . 

In fact, the trend of increasing differential abundance 
with increasing condensation temperature approximately 
mirrors th e pattern of depletion of elements in ISM gas 
(jYinll200^ . This simply reflects the preferential conden¬ 
sation of refractory elements onto interstellar grains and 
the universality of chemistry. If interstellar dust does not 
perfectly follow the gas, this will produce variation in el¬ 
ement abundances that are proportional to the extent of 
gas-dust segregation and the dust composition. All else 
being equal, the abundances of the more refractory ele¬ 
ments (which are mostly in the dust) would be expected 
to vary more than volatile elements (which are mostly in 
the gas), in a manner that roughly correlates with Tq. 

There is evidence for, and theo retical predict i ons of , 
dust-gas segregation in the ISM. iPadoan et al.l (|2006f l 
found that the power-law describing the spatial distribu¬ 
tion of extinction (a tracer of the column density of dust) 
on the sky is shallower than that of (a tracer of the 

gas column density), such that in the densest regions of 
the ISM, extinction by dust is less than that predicted 
by perfect correspondence with the gas. This could be 
explained by grain growth in excess of that predicted by 
models. (CO c ondensation would p resum ably act in the 
oppos ite sense IWhittet et al.l 1201011 . But IPadoan et ahl 
(|200(ill offer spatial variation in dust-to-gas ratio as an 
alternative explanation. 

Dust-gas segregation could occur at three different 
scales; (i) over the extent of a giant molecular cloud or 
star-forming region; (ii) in cloud cores that collapse to 
form individual stellar systems; and (iii) in the acc retion 
disks around young stellar objects. IPraind (|2011l l esti¬ 
mated that dust drift due to radiation pressure from O 
stars can remove dust from the centers of H II regions 
in < 1 Myr, but since such low-density regions are not 
themselves the site of star formation, this effect should 
not manifest itself in the relative abundances of stellar 
photospheres. 


iBellanI (|2008ll showed that the different dynamics of 
dust and gas during Bondi-type accretion flow can en¬ 
hance dust-to-gas ratios by an order of magnitude in 
cloud cores. The mechanism modeled by IBellanI (120081 1 
relies on dust velocities of several km sec“^ relative to 
the gas around the growing cloud core. The size distri¬ 
bution of interstellar dust grains peaks sharply near 0.2- 
0.3 fj,ia (|Drainell200,3jl . These grains may be dynamically 
decoupled from the diffuse interstellar medium (number 
density n ^ O.I cm“^, stopping distance ^ 40 pc), and 
there is tentative evid ence for this in th e trajectories 
of interstellar meteors (iTavlor et al.llIOOli ) and Ulysses 
spacecraft measurements (iKriiger fc Grunll20d^ . How¬ 
ever, these grains will be tightly coupled to the molecular 
gas (n ^ 10^ — 10^ cm“^) surrounding a cloud core, i.e. 
over length scales <§; 1 pc and much smaller than a typical 
cloud size. Under such conditio ns it seems unlikely that 
dust acceleration mechanisms (|Yanl 120091: iHoang et ahl 
IMl e.g.,) can achieve equipartion of energy between 


ISM and dust and speeds of ^10 km s“^. 

The second scenario which could produce variation in 
the dust-to-gas r atio in cloud cores is drift in duced by ra¬ 
diation pressure. iWhitworth &: Bat^ (|2002f) showed that 
the inward radial drift of 0.1-/rm grains under the influ¬ 
ence of a typical radiation field can treble the dust con¬ 
centration in a static gas sphere of a few solar masses at 
10 K in 10 Myr. The effect scales linearly with the inten¬ 
sity of the external radiation field and the inverse of the 
characteristic cloud column density. Grains reach termi¬ 
nal velocity on a timescale of 10^ yr, i.e. much shorter 
than the cloud lifetime. Thus any density enhancement 
can be expected to grow linearly and the total dust en¬ 
hancement can be expected to scale with the total ex¬ 
ternal r adiation experien ced by the cloud core over its 
lifetime. iSeo et aH (1201111 numerically simulated this ef¬ 
fect, including gas dynamics as well as coupling between 
the gas and dust, and found that the dust concentration 
is enhanced by about an order of magnitude in a narrow, 
inward-propagating shell, and depleted exterior to that 
shell. While the mean dust-to-gas ratio of a cloud is not 
changed by migration, truncation of the cloud by photo¬ 
evaporation or internal collapse (to the exclusion of outer 
regions) would produce a me tal-enhanc e d obj ect. Both 
IWhitworth fc Bat^ (j2002l l and lSeo et al.l (|20111 ) point out 
the relevance of this process to the metallicities of stars 
and the formation of their planets. 

A third scenario for dust-gas segregation involves 
accreting protostars and their disks, or concomi- 
tant Bondi-Hoyle accr etion from the molecular cloud 
(|ThrooD &: Ballvl 1200811 . Photoevaporation of gas, but 
not dust, from a disk produces an enhancement in the 
dust-to-gas ratio which is inherited by the star as disk 
accretion continues. Photoevaporation of disks driven 
by X-rays from the central star has been proposed to ex- 
plain the final, rap id stage of circumstellar disk clearing 
(lOwen et al.ll201 2l). Observatio ns of blue -shifted lines of 
Ne II ( Pascucci fc Sterzi3l2009ll and O I (jHartigan et ahl 
Il995f) indicate that heavy volatile elements as well as H 
and He are lost in these winds. The temperature of the 
X-ray-heated “surface” of the disk fr om which the wind s 
flow is heated to a few thousand K (lOwen et al.ll20I^ . 
but gravitational settling (and the formation of planetes- 
imals) keeps dust near the cooler mid-plane. Disk masses 
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are typically 1% of th e central star, with a large scatter 
(|Andrews et al.l[2013[) : assuming the convection zone of 
a young solar-mass star contains O.IMq, accretion of an 
entire disk which has been severaly gas-depleted would 
increase relative elemental abundances by < 0.04 dex, or 
< 0.01 dex if the convection zone contains O.4M0. This 
may fall short of explaining some of the observations. 

Two important clues to the mechanism(s) of the ob¬ 
served variation are the spatial scale on which the vari¬ 
ation occurs, and the Tq of the “elbow” in the rela¬ 
tive abundance variation below which abundances do 
not vary in any systematic manner. Surveys of nearby, 
solar-age field stars offer little information on the spa¬ 
tial scale of the segregation as the stars are far from 
their birthplace. lOnehag et al.l (|2014D found no signif¬ 
icant difference between M67 cluster stars and the Sun, 
and proposed that the Sun also formed in a dense clus¬ 
ter analogous to M67. They suggested that the common 
pattern of relative abundances was established by the re¬ 
moval of dust by radiation drift. But there are several 
problems with this explanation. Star formation does not 
take place in HII regions but in surrounding neutral gas, 
and, as they point out, the expansion of an HII region 
should outpace dust drift. This explanation also requires 
most stars to form after an earlier generation of massive 
stars that move the dust. It also requires that the vast 
majority of nearby solar-type stars have not formed in 
dense stell ar clusters, something not supported by cluster 
statistics (|Williams fc Gaidosll2007ll . Finally, abundance 
differences between stars in physical pairs requires gas- 
dust segregation on scales smaller than the cloud core. 

The other clue is the value Tq below which relative 
abundance variation disappears, as this suggests the tem¬ 
perature and hence location where the dust and gas seg¬ 
regation occurs. If this “knee” in the abundance vari- 
ation vs. Tg is 100 0 K, as many data sets suggest 
(|Melendez et al.1 [200911 . this would seem to rule out a 
molecular cloud setting with T ~ 20 — 100 K where even 
some G and O condense as organics and ices. All of these 
observations point to gas-dust segregation during the for¬ 
mation or subsequent evolution of an accretion disk as a 
plausible explanation of the abundance trends. 

6. DISCUSSION 

Neither observations of M dwarfs nor models of Galac¬ 
tic chemical evolution support the premise that there 
are a significant number of stellar systems with primor¬ 
dial C/0 ~ 1 in the solar neighborhood, i.e. within the 
samples of Doppler radial velocity surveys. Moreover, if 
the process of planet formation is universal, then studies 
of primitive meteorites show that the chemical composi¬ 
tion of small, rocky planets is controlled largely by the 
composition and thermal processing of pre-existing dust 
from the parent molecular cloud, rather than condensa¬ 
tion from gas with the stellar composition, and that this 
dust is likely to have a C/0 reflecting that of the bulk 
ISM. Finally, small differences in the relative patterns be¬ 
tween the Sun and solar “twins” can best be explained 
by dust-gas segregation at stellar scales and tempera¬ 
tures of up to 1000 K, i.e. in accreting protostars, rather 
than planet formation per se. 

When comparing the Sun to the solar neighborhood 
it is usually assumed that they share a common his¬ 
tory. However, it is also possible that the two are 


not related. iNieva fc Przvbili^ (I2012D estimated “cos¬ 
mic” standard abundances (rather the present-day abun¬ 
dances in the solar neighborhood) using bright, slowly- 
rotating, early-type B stars. Due to their short main se¬ 
quence lives, B stars contain abundance information that 
is both contemparaneous and local, and their purely ra¬ 
diative atmospheres are comparatively simple to model, 
although UV photoionizati on rates must be correctly 
modeled (iLvubimkovI l2013l l. INieva fc Przvbill^ (1201^ 
report metallicities that are very close to solar, depsite 
4.5 Gyr of intervening GCE, and a mean C/0 of 0.37, 
lower than the solar value of 0.55. They explain both 
discrepancies by appealing to outward migration of the 
Sun in the galactic disk. These findings warrant further 
investigation. 

Depletion patterns in the ISM suggest that the C/0 
of dust largely reflects the ISM; O is incorporated more 
rapidly in grains than C in molecular clouds but it 
is depleted from grains more quickly in the intercloud 
medium, and these two effects approximately balance. 
Variation in the efficiency of grain erosion by SN shocks 
and UV radiation could produce modest variations in 
the C/0, however. While the phases of O in interstel¬ 
lar grains are clearly water ice and silicates, the phases 
of C are more controversial. Volatility patterns suggest 
C is not in refractory phase(s) like graphite and carbides 
but instead more volatile organic molecules, in particular 
aromatic hydroca rbons (j.Ione i 120091 : l.lones et 'm 120131 : 
iGhiar et al.l l2013f) . Aliphatic hydrocarbons have been 
suggested as the source of the ubiquitious ISM absorbtion 
feature at 3.4/im but the absence of polarization in this 
line is difficult to reconcile with a scenario of c ondensa¬ 
tion o f carbonaceous mantles onto silicate cores (iLi et al.l 
l2014|l . Surface reactions could explain why growth rate 
of G in grains is very slow with respect to other elements. 

While these considerations plus a bulk ISM C/0< 1 
are sufficient to explain the oxidized nature of dust in 
the Solar System, it remains to be explained why prim¬ 
itive meteorites in the Solar System are very depleted 
in both G and O, especially G. Based on the mete- 
ori te abundances and associated uncertainties compiled 
bv iLodd^ (120031) and solar abundances and uncertain- 
ties from iCaffau et al.l (120111) . the depletion has a sto- 
chiometric ratio of 0.93 ±0.34, i.e. consistent with unity 
and depletion by the formation and removal of GO. In 
this scenario the removal of G or O is not controlled by 
volatility, otherwise either G or O would be far more 
depleted than the other, depending on the effective Tq- 
This can occur if the reaction occurs in situ, i.e. at tem¬ 
peratures lower than the Tcs of the G and O phases 
(~ lOOK) where both elements are retained in grains, 
and/or where vapor is not lost by turbulent mixing, i.e. 
in “dead zone” where turbulence is suppressed. 

Primitive meteorites, considered analogs to the now- 
lost building blocks of the Solar System, also include 
the highly-reduced enstatite chondrites. Enstatite chon¬ 
drites contain carbide minerals but virtually no water 
and have been explained b y equilibration of solids with 
a gas having C/0 = 0.83 (|Grossman et al.l[2008[) . In¬ 
terestingly, removal of an amount of C and O equal to 
that incorporated in Cl chondrites produce s a re sidual 
gas with a C/0 of 0.93. iHutson fc Ruzick^ (|2000D have 
proposed that removal of refractories from a gas of solar 
composition and equilibration of solids with the remain- 
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ing gas could produce the precursors for enstatite chon¬ 
drites. If enstatite chondrites formed from recycling of 
the precursors of terrestrial planets in an oxygen-poor gas 
it could explain why the enstatite chondrites and Earth 
lie along the same mass-dependent fractionation line in a 
three oxygen-isotope plot. This equilibration would nec¬ 
essarily have occurred at a comparatively low dust-to-gas 
ratio and lower temperatures where moderately volatile 
elements such as Fe would not be lost. 

Al-Mg isotopes of CAIs in unequilibrated enstatite 
chondrites indicate that these refractory inclusions 
formed in the same region as CAIs in other chondrites, 
but were subsequ ently exposed to reducing conditions 
(|Guan et al.ll2000f) . Reducing conditions could have been 
established sequentially or in a different part of the pro¬ 
toplanetary disk. The formation times of enstatite chon¬ 
drites is an area of active research. Dating of sulfides in 
unequilibrated E chondrites using the Fe-Ni and Mn-Cr 
short-lived radionuclide chro nometers gives ages of 12- 
13 Myr after CAI formation (iWadhwa et al.l[l997fl . Re¬ 
analysis of the Mn-Cr isotope data for a sin gle sulfide in 
the MacAlpine Hills 88136 EL3 chondrite ( Guan et al.l 
l2007t I Telus et all l2012t) pl us an initial ®^Mn/®®Mn of 
5.1 X 10“® (|Yin et aP 1200711 gives an age of 10 Myr. 

In summary, I propose the following scenario for the 
chemistry of rocky planets: 

• Planet formation proceeds from a two-component 
ISM (gas and dust) which are never fully equili- 
braxted during the formation process. 

• The composition of the precursor material of Solar 
System material is interstellar dust with refractory 
abundances and a C/0 ratio approximately equal 
to the bulk ISM; this was set by grain growth in 
molecular clouds and destruction between clouds. 

• Planetesimals, as represented by carbonaceous (Cl) 
chondrites, formed under the oxidizing, high dust- 
to-gas ratio conditions established by settling of in¬ 
terstellar grains to the disk mid-plane and subse¬ 
quent depletion of C and O by formation and re¬ 
moval of CO. 


• Enstatite chondrites represent reduced material 
that equilibrated with gas after removal of solids 
in approximately Cl proportions, possibly during a 
second generation of planetesimal formation. 

• The chemical composition of rocky exoplanets 
could be set by mixing of oxidized and reduced gen¬ 
erations of planetesimals. The mixing ratio could 
be determined by the efhency with which the first, 
oxidizing generation incorporated disk solids, as 
well as the dynamics of accretion in the disk. 


Giant planets, unlike the small, rocky planets consid¬ 
ered in this work, accrete massive gas envelopes that 
would include volatile species such as CO, and the con¬ 
siderations described above do not preclude the possi¬ 
bility of C-rich atmospheres in giant planets. Detec- 


versial (e.e. iMadhusudhan et al.l 

20111 iCrossfield et al. 

120121 ISwain et al.ll2013l:lLine et al. 

12014 IStevenson et al. 


12014 iHansen et al.l 12014) . If the core-first model of gi¬ 
ant planet formation is correct, removal of O as silicates 
and sequestration into a core increases the C/0 ratio of 
the gas that is subsequently captured into the planet’s 
envelope. Beyond the ice-line, condensation of water and 
removal of more O would drive the C/0 of the gas even 
closer to unity (|Oberg et al.ll20lil) . Dissociation of CO 
in the disk and removal of the O as water ice would fur¬ 
ther enhance this ratio, and indeed the enrichment of 
planetary water in ^^O and ^®0 relativ e to the Sun i s 
thought to be a signature of this process l|Clavtonll2002D . 
These processes may facilitate the formation of reduced, 
enstatite chondrite-like planetesimals as discussed above. 
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APPENDIX 


GALACTIC CHEMICAL EVOLUTION MODEL 

I computed changes in the abundance of five isotopes in the vicinity of the solar galactocetnric radius: 

^®0, and ^®0. The model accounts for the production of these isotopes and their release into the ISM in SN 
explosions, winds from massive stars, and AGB winds, incorporation of isotopes into long-lived low-mass stars, and 
dilution of the ISM by the infall of metal-poor gas. The ISM is described by a two-box model with a lower-density, 
warmer intercloud medium (ICM) which spawns molecular clouds, and a giant molecular cloud (GMC) component 
that can form stars. The mean lifetime of molecular clouds is tgmc and the residence time of gas in the IGM before 
condensing into a cloud is tjcm • The equations of motion for the mass surface density m of the two components are: 


dmicM micM , niGMC 


dt 


"TICM 


tgmc 


+ ‘S'icM + ‘S'agb + F, 


drriGMC micM »71gmc 


dt 


TICM 


tgmc 


+ -ScMC — RM, 


(Al) 
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where F is the infall rate, R is the rate of formation rate of stars with M > IMq, A4 is the ratio of the total stellar 
mass to the mass in stars > IMq, S'icm is the flux of SN ejecta into the IGM from progenitors that explode in time 
T > tgmc, S'gmc is the flux of SN ejecta into the parent molecular cloud from progenitors that contribute to GMGs 
in time T < tgmc, and S'agb is the wind from AGB stars, all contributing to the ICM. The fluxes from AGB winds 
and SN ejecta are given by: 

^ (A3) 

J Ml ) 

where (M) is the mean mass of stars with M* > IM 0 , E is the ejected mass in winds and/or explosions, / is 
the fractional number of stars per unit mass in the IMF, and Mi and M 2 are the minimum and maximum masses 
contributing mass to the three different cases, and P is the mass in the wind or ejecta from a progenitor of mass M*. 
The equations for the mass surface densities of stars (m*) and stellar remnants {uir) are 
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dt 
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where the limits of integation are over any stars that are moving off the main sequence. The star formation rate is 
related to the total mass surface density of gas using a Schmidt-Kennicut law: 
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The equations governing the the mass fraction of the ith isotope is: 
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Aq are the isotopic abundances of the infalling gas, and 
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Progenitor mass (solar mass) 

Figure 10. Predicted mass of oxygen (blue, solid) and carbon (red, dashed) from stars vs. progenitor mass for a solar-metallicity 
population weighted by an IMF with power-law index a. = 2.6 

where P® is the production (equal to the nucleosynthetic yield plus the original mass of the isotope). The C/0 ratio 
is calculated by summing the appropriate mass fractions divided by the atomic weights. An approximate metallicity 
is also calcated from the sum of all the isotopes relative to the abundances in the Sun. 

I considered nucleosynthesis in the p rogenitors of AG B stars (1-6. SMq), “super-AGB” stars (6.5-11Mq), and SN (11- 
I 2 OM 0 ). 1 used the AGB yields from iKarakasI (120101) for AGB progenitors up t o 6.5M f;). For super-A GB progenitors 
with masses between 7 and llM f^ 1 used the production s fro mlPohertv et all (I2014D or iSiessI (I2010D otherwise. For 
SN yields 1 adopted the values in iKobavashi et ahl (l2006f ) and iKobavashi et al.l (I2011D for regular SN with energies of 
10®^ erg and their re-run calcalation s for the cases o fISM q and Z = 0.004, and M* = 25 and Z = 0.02. These 
were supplemented with values from iPortinari et al.l (Il998f ) for progenitor masses of 12, 60, 100, and 120Mq. In order 
to arrive at a solar C/0 with reasonable cho ices of parameters 1 fin d it necessary to include yields from the winds of 
massive stars with M* > 4 OM 0 estimated bv iPortinari et ahl (jl998D (Table 3 in that work). 

The literature for SN and AGB nucleosynthesis is very heterogeneous, with calculations performed for different 
ranges/values of progenitor masses and metallicities, and it was necessary to interpolate or extrapolate for some 
val ues. Spec i fically , to estimate the production of 6 . 5 M 0 AGB stars with sub-solar progenitor metallicities 1 scaled 
the IKarakasI (120101) values for sub-solar metallicity 6 M 0 progenitors by the ratio of the 6 . 5 M 0 to 6 M 0 production for 
solar-metallicity. For productions from super-AGB stars with progenitor masses outside the caclulated range, 1 either 
used the production from the most massive progenitor with the same metallicity, or linearly interpolated between the 
values for the most massive AGB and least massive super-AGB progenitors with the same metall i city. F or a review 
of the many parameters and uncertainties that enter these calculations see iKarakas fc Lattanziol (I2014D . 1 linearly 
interpolated the production onto a grid of 1000 masses over I-I 2 OM 0 with intervals chosen such that the IMF has 
equal total mass in each bin. Figure [10] shows the calcu late d yi elds (production - initial incorporation) for solar 
metallicity. These values were then used to evaluate Eqns. IA3l and lA9l 

1 searched for the parameter values that best reproduce the observational constraints using Monte-Garlo Markov 
Ghain analysis. The constraints are the age of the Sun, the present mass surface densities of stars, stellar remnants, 
and total gas at the solar galactocentric radius, solar C/0, the current metallicity of the ISM (taken to be -1-0.14 
based on the metallicities of M44 and the Hyades), and t he mean (-0.05) and intrinsic standard-deviation (0.18 dex) 
of the metallicity distribution in the solar neighborhood ([Gaidos et al.ll201^ . 1 adopted standard errors of 0.12 in 
C/0 (jCaffau et al.ll20lil) . 1.5, 1, and 1 M 0 pc“^ for stars, remnants, and gas, 50 Myr for the absolute age of the Sun, 
and 0.01 and 0.02 dex in the mean and standard deviation of the present metallicity. 1 varied the exponential infall 
timescale, the metallicity of the infalling material, the age of the Galactic disk in the solar annulus, the high-mass 
IMF index a, and the Schmidt-Kennicut index (3. 

The simulation is able to adequately reproduce the observed properties of the Galactic disk in the solar neighborhood 
(y^ = 35 with v = i). The most significant deviation is the over-prediction of the present ISM gas mass and 
underprediction of the current metallicity. The adopted or fit values for the simulation are given in Table [I] The 
predicted evolution of the SFR, mass surface densities, and the metallicity of the ISM are pl otted in Fig. ITT] Estimate s 
of the current total star formation rate of the Milky Way cluster around 1.9 ±0.4 M 0 yr“^ (iChomiu k &: Povichl 20111) . 
Presuming that star formation has an exponential radial distribution with scale length of 3 5 knc ([Wolfire et al. I20(i^ 
then the rate at the solar radius will be about 2.7 ± 0.6 x 10“^ M 0 Myr“^ pc“^. This compares favorably with the 
predicted value at the present (Fig. ITT]) . 
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Figure 11. Predictions of the GCE model compared to the Sun and values at present. Top: star formation rate per area compared to 
the present estimated value. Middle: mass surfa ce densities of stars (sol id, red), gas (green, dotted), and stellar remnants (blue, dashed) 
compared to estimates of present values. fDoints. INaab &: Qstrik^l20Q6i') . Bottom: metallicity, here taken to be [C+O/H]. 


Table 1 

GCE Simulation Parameters 


Symbol 

Parameter 

Value 

Units 

Reference 

Solar Parameters 


Solar oxygen 

6.73 X 10-^ 

— 

Caffau et al. f20111 

XC 

Solar carbon 

2.73 X 10-3 

— 

Caffau et al. f20lG 


Stellar Parameters 



a 

IMF index 

2.61 

— 

Fit 

M 

Low mass/high-mass ratio 

2.00 

— 

Fit 


Galactic Parameters 



m,(T) 

Present stellar mass density 

28.5 ± 1.5 

M0 pc-2 

Naab & Ostriker ('2006) 

mg(T) 

Present gas mass density 

13.5 ± 1 

Mq pc-3 

Naab Sz Ostriker ('2006') 

mr{T) 

Present remnant mass density 

6±1 

M© pc“^ 

Naab & Ostriker ('2006') 

'TGMC 

GMC gas residence time 

30 

Myr 


'QMC 

IMG gas residence time 

150 

Myr 


Tdisk 

Age of disk 

9554 

Myr 

Fit 

"^infall 

Infall e-folding time 

2630 

Myr 

Fit 

[M/H]o 

Metallicity of infalling gas 

-1.3 

— 

Fit 

P 

SFR rate index 

1.71 

— 

Fit 



















